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Abstract 

> 

Recent  backscatter  results  generated  by  Tappert,  using  parabolic  equation 
and  Monte  Carlo  methods,  afford  an  excellent  opportunity  to  assess  the  ade¬ 
quacy  of  concepts  such  as  shadowing  and  surface  tilting  often  used  in  concert 
with  ray  concepts  for  line  of  sight  backscatter  calculations.  In  this  study, 
results  of  first  order  scatter  from  rough  surfaces  are  used  in  conjunction 


with  waveguide  formalism  to  calculate  clutter  from  distant  ranges  in  tropos¬ 
pheric  waveguide  environments.  Comparisons  are  made  with  Tappert’s  results  at 
9.6  GHz  for  the  standard  atmosphere,  and  for  14  and  28  m  evaporation  ducts. 
Results  apply  to  wind  speeds  of  10,  20,  30  and  40  knots.  Averaged  backscat- 


tered  signals  calculated  by  the  two  methods,  for  a  transmitter  altitude  of  25 


m,  agree  to  within  about(d£  10  dB.  This  is  considered  surprisingly  good  agree¬ 
ment  in  view  of  the  many  uncertainties  and  approximations  involved  in  the 


1.0  INTRODUCTION 


'ff 

Tappert  ,  in  a  recent  study,  has  used  parabolic  equation  and  Monte  Carlo 
methods  to  calculate  sea  clutter  for  a  variety  of  sea  states  and  evaporation 
ducts.  In  particular,  he  has  carried  the  clutter  calculations  well  beyond  the 
normal  horizon.  An  interesting  question  is  whether  waveguide  methods  might 
also  be  useful  for  such  calculations?  An  immediate  response  would  quite  pos¬ 
sibly  be  no,  since  waveguide  modes  are  generally  complex  and  there  would 
seemingly  be  no  way  of  accommodating  such  quantities  within  the  conventional 
framework  of  line  of  sight  rough  surface  scatter  which  deals  with  real  angles 
of  incidence.  Nevertheless,  Tappert 's  work  has  motivated  the  present  ex¬ 
ploratory  study  of  sea  clutter  via  normal  mode  methods. 

In  the  following  section,  abundant  useage  is  made  of  first  order  scatter 
theory  from  a  rough  surface  as  developed  by  Ulaby,  Moore  and  Fung^ .  In  par¬ 
ticular,  first  order  theory  is  used  in  conjunction  with  normal  mode  formalism. 
As  will  be  shown,  this  leads  to  the  physically  satisfying  result  of  the  back- 
scattered  power  falling  off  inversely  with  range  along  with  an  attenuation 
associated  with  two  way  transit  between  transmitter  and  receiver.  First  order 
scatter  theory  is  applied  to  capillary  waves  on  an  otherwise  flat  surface.  In 
actuality,  the  capillary  waves  ride  on  longer  wavelength  waves  which  scatter 
energy  out  of  the  specular  direction  and  thereby  augment  modal  attenuation 
rates.  Approximate  allowance  for  this  effect  is  made  by  assuming  the  surface 

obeys  Gaussian  statistics  and  by  using  Kirchhoff  theory  to  modify  the  smooth 

2  3 

surface  ground  reflection  coeffi  lai.r.  ’  Additionally,  the  long  wavelength 
waves  effect  a  tilting  of  the  surf..  ,  thereby  altering  the  local  angle  of  in¬ 
cidence.  Moreover,  the  longwaves  shadow  portions  of  the  scattering  surface. 
Awareness  of  these  effects  has  given  rise  to  the  so  called  "composite  surface 
model"  which  has  been  used  extensively  for  within  the  line  of  sight  calcula¬ 
tions  in  non-ducting  environments .  ^ Though  lacking  in  rigor,  these  long 


2 

*  This  citation  refers  to  a  hand-written  report  prepared  by  Dr.  F.  D.  Tappert 
of  the  University  of  Miami,  while  working  at  NOSC  during  the  summer  of  1989. 

For  further  information,  please  contact  R.  A.  Pappert,  Code  542,  or  H.  V.  Hitney, 
Code  543,  Naval  Ocean  Systems  Center,  San  Diego,  CA. 


wavelength  effects  are  treated  in  section  III  by  the  same  methods  convention¬ 
ally  used  for  line  of  sight  calculations  in  non-ducting  environments. 

Comparisons  with  Tappert's  results  are  given  in  section  IV  and 
a  general  discussion  in  section  V  concludes  the  paper. 

2.0  BRAGG  SCATTER 

In  the  following  development,  a  rectangular  coordinate  system  is  used 
with  x  the  range  variable,  z  the  altitude  variable  with  z  —  0  corresponding  to 
ground  level.  A  thin  vacuum  layer  is  assumed  immediately  above  a  random 
rippled,  but  otherwise  flat,  ocean  surface  and  the  polarization  is  taken  to  be 
horizontal.  Ulably  et  al . ^  have  treated  in  considerable  detail  the  problem  of 
first  order  scatter  from  such  a  surface  for  plane  wave  incidence.  Their 
result  for  the  spectral  decomposition  of  the  scattered  field  due  to  a  unit 
amplitude  incident  plane  wave  in  the  x-z  plane  is  (a  time  dependence  exp(jwt) 
is  assumed) 

EyS)  ~  fc  X/Uy ( kx , ky ) exP ( J kxx+ j kyy - j kz z ) dkxdky ,  (1) 


where 


U(kX’V  “  -h(kX>ky)f. 


h<kx-y  -  <kzkz+kx>/D- 


kz  -  <k2-kH)1/2.  <  -  (k2n2-k2-k2)V2, 


d  -  (vyoc^v. 


(2) 

O) 

(4) 

(5) 
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(6) 


f  -  jk2(n2-l)(l+Ri)z(kx+kS,ky), 

2  a _  e_ 

g  i^£0  V 

C  -  (rig-S2)1/2 

1  °° 

Z(kX,ky)  -  2^//z(x.y)exp(-jkx-jkyy)dxdy. 

S  -  sin# ,  C  -  cosd , 


0  -  plane  wave  angle  of  incidence, 


k  -  free  space  wave  number, 


a)  -  circular  frequency, 


a  -  ground  conductivity, 


—  -  ground  dielectric  constant, 
£0 


z(x,y)  -  random  function  defining  rippled  surface. 


(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

(14) 

(15) 


Unless  otherwise  specified,  MKS  units  are  assumed.  In  the  above  Rj^is  the 
Fresnel  reflection  coefficient  for  horizontal  polarization.  Signs  of  the 
square  roots  which  occur  in  the  above  equations  are  taken  in  such  a  way  as  to 
assure  either  upward  propagation  or  evanescence  for  z>o  and  downward  propaga¬ 
tion  for  z<o. 
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The  broadside  electric  field,  in  the  assumed  vacuum  region  bordering  the 
surface,  generated  by  a  horizontal  electric  dipole,  p  ,  may  be  written  as  the 
following  modal  decomposition  (kx»l  is  assumed) 


'y110  “J  gVym(zT)exP(-jkSmx)(exP(JkCraz)+  RimexP( ' jkCmz) >  * 


(16) 


where  m  is  a  modal  index,  e 

and  is  normalized  to  1  +  Rlm 
given  by 


(z~)  is  the  modal  height  gain  at  the  transmitter 
ym  i 

at  the  ground.  The  excitation  factor,  A  ,  is 


(17) 


where  the  modal  function,  F,  is 


F  -  1  -  RjRj.,  (18) 

and  R  is  the  plane  wave  reflection  coefficient  for  horizontal  polarization 
from  everything  above  z  -  0  with  vacuum  below  it.  The  constant  A  is  given  by 


A  “  4 ^  tk)1/2exP(J*/4).  (19) 

o 

To  good  approximation,  the  factor  1+Rj^  which  appears  in  Eq.(6)  may  be  written 
as 


1  +  Rj_ 


2C 


^1/r 


(20) 
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The  amplitude  at  the  surface,  of  the  incident  plane  wave  associated  with  the 
dipole  generated  field,  is  interpreted  to  be  the  factor,  apart  from  the  range 
dependent  exponential  term,  which  multiplies  the  exp(jkCmz)  term  in  Eq.(16). 
Combining  this  factor  with  Eq.(l)  and  using  Eqs.(9)  and  (20)  yields 


,(s)=  _  jAkW-l)1/2 


&_ _  JTJjMkx,k  )z(x  ,y  )exp(jkx(x-x  )+jk  (y-y  )-jkz(z-z  )) 


2n 


(i01/2 


C  A  exp(-jkSx  )e.(zT)  dkdkdx  dy 
m  m  r'  J  m  '  ym  T'  x  y  J 


(21) 


To  arrive  at  Eq.(21)  it  has  been  tacitly  assumed  that  the  incident  radiation 

t 

is  an  outward  propagating  pulse  with  x  much  greater  than  the  pulse  length  so 

that  the  geometrical  factor  1/Vx7  has  been  taken  to  be  constant  over  the  il¬ 
luminated  area. 

The  integrals  over  k  and  k  are  now  transformed  to  cylindrical  coor- 

x  y 

dinates  via  the  transformations 


kx  -  ksin0cos<£,  k^  -  ksin^sin^, 


kz  -  kcostf,  dk^dky  -  k  sin0cos0d<£ , 


:-x  -  rcos^r,  y-y  -  rsin^r> 


Also,  to  good  approximation 


Mk„  k  )  -  4-  — 1 


*■  v  k'  k^2-!)1/2- 

z  v  g 


(22) 

(23) 

(24) 


(25) 
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so  that  the  scattered  field  becomes 


3  oo  '  '  2+jc0 

EyS)“-jf2  fS  dx'dy'^^SCmAmexP('jkSmx'>eym(ZT)  JoSCexP< ‘ JkCz)d* 


y  2^ 


•  J^expCjkrScos  )  )  d<f> . 


The  upper  limit  of  the  theta  integration  has  been  determined  by  the  condition 
kz~>  -j«>  as  kx,k  -+  00 .  The  integral  over  <f>  is  proportional  to  the  zero  order 
Bessel  function,  J  ,  so  that  the  scattered  field  becomes 


3  *  .  . 


J.fdx'd y 


z(x' , y ' ) 


f+jco 


y  n 


SCmAmexp(‘jkSmx  )eym<zT)/  SCJQ(krS)exp( - jkCz)d0 . 


The  integral  over  9  can  be  written  in  terms  of  an  integral  over  the  Hankel 
function  of  the  second  kind,  .  Effecting  that  transformation  allows  for 


the  scattered  field  to  be  written  as 


f+j~ 


E<S)  =  -jf£  jjdx'dy'  _  SCmA^exp(  r  '„x  '  )cy.m(  zT)f  SCH<2 }  (krS)exp(  -  jkCz)d0 . 
--  -w  jx'  m  n  . 

.-.JCO 


Equation  (28)  represents  the  scattered  field  in  tb«  ’-bin  vacuum  region  in 
the  neighborhood  of  the  surface  without  full  allowance  for  stratification. 
This  field  is  generalized  to  fully  apply  to  a  stratified  media  by  allowing  for 
the  repeated  reflections  which  occur  between  ground  and  the  stratified  media. 
To  make  this  more  precise,  a  vacuum  will  be  assumed  between  ground  level  and 
an  infinitesimal  level  a.  The  plane  wave  reflection  coefficient  from  the 
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ground  is  simply  the  Rj^  coefficient  already  introduced.  The  plane  wave 
reflection  coefficient  from  everything  above  level  a  and  referenced  to  level  a 

will  be  denoted  by  R^q.  The  infinite  sum  of  reflections  associated  with  the 
process  whereby  the  upgoing  wave  from  the  surface  goes  directly  to  the 
receiver  at  z<a  or  first  reflects  off  the  upper  stratification  and  then  off 
the  ground  to  the  receiver  etc.  leads  to  the  expression 

z  -  a  _  _ 

z  +  z  +  --  ->exp(  -  jkCz)/(l-Rj^  R^Qexp( -2jkCa) )  .  (29) 

z  -  o  _  _ 


The  infinite  sum  of  reflections  associated  with  the  process  whereby  the  upgo¬ 
ing  wave  first  reflects  from  the  stratification  to  the  receiver  and  then  from 
the  stratification  to  the  ground  and  once  again  from  the  stratification  to  the 
receiver  etc.  yields 


z  -  a 


R  Qexp(-jkC(2a-z))/(l-R_L  Rj.aexp( -2jkCa) ) .  (30) 

When  the  z  dependent  term  in  the  theta  integral  of  Eq.(28)  is  replaced  by  the 
sum  of  Eqs .  (29)  and  (30),  the  limit  a,  z-*0  taken  and  the  theta  integral 
evaluated  by  residue  theory  there  results 


E<s>-  Ak3//  dx'dy,^^5C«V’‘P(JkSm’‘'>V<ZT)SSnCnVo2)<kV)<1+Rlm)-  <31> 
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The  ter.,-  x  +  R^n  is  simply  the  modal  height  gain  at  the  ground.  For  backscat- 
ter  at  the  transmitter  location,  the  1  +  R^n  term  need  only  to  be  replaced  by 

the  modal  height  gain  e  n(zT)  with  the  understanding  that  e  (0)  -  1  +  R^n- 

*  (2) 
Making  this  replacement  as  well  as  the  asymptotic  replacement  for  and 

noting  that  Sn  <*  1,  Eq.(31)  becomes 


{t.\  ®  .  ,  z  (x  ,y  )  ,  i 

y  )=  dx  dy  / — —  §CmVxp('jkSmx  )eym(zT)  r  £CnAnexp(  '■ jkSnr)eyn(zT)  • 

y  x  'r 


where 


B  -  ^  2k5. 
o 


It  can  further  be  shown  that  if  kA<£y«l  where  is  a  measure  of  the  angular 

beam  pattern  (assumed  to  be  on  the  order  of  hundreths  of  a  radian)  and  y  a 
measure  of  the  correlation  length  of  the  capillary  waves  (presumed  to  be  on 

I 

the  order  of  centimeters)  then  to  good  approximation  r  can  be  replaced  by  x 
in  Eq . ( 32 )  so  that 


EyS  =*  BJ7  dx'dy' 
y  -co 


FmAmexp(-jkSmx  )eymUT)  £CnVxp( ‘ JkSnx  )eyn(zT)  •  (34) 


Equation  (34)  is  the  y  component  of  the  electric  field  at  the  receiver. 
To  calculate  the  Poynting  flux  the  Hz  field  is  also  required.  From  Maxwell's 

curl2  equation  the  Hz  field  is 


-  E<SV 
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where  tj  is  the  free  space  impedance,  JnQ /e  ,  with  nQ  the  freespace  perme¬ 
ability  and  £q  the  free  space  permittivity.  The  time  averaged  Poynting  flux 
is  (the  *  denotes  the  complex  conjugate  and  the  braces,  <> ,  an  ensemble 
average) . 


2  <  E<S)  H<s>*> 


ml 

2r7 


r  r  f  9  9  9  99  9  9  9  9  99  _ = _ 

JJJJ  dx  dy  dx  dy  <z(x  ,y  )z(x  ,y  )>  ,  , , 

‘00  XX 


C  C  C  AA  A  A  exp (jk(  (S  +S  )x  -(S  +SJx  )) 
n*ripq  ranpqmnpq  rvj  n  nr  N  p  q'  '' 


eym(zT)eyn(zT)eyp(zT^eyq^zT) ' 


(36) 


The  space  integrals  in  Eq.(36)  are  next  transformed  as  follows 


_  II  1  _  t  1  1 

X_X  -  X  ,  y  -  y  -  y, 

1  III  1  ill 

X  =  2  (x  +x  )  ,  Y  -  (y  +y  )  , 

(37) 

1  1  II  It  _  _ 

dx  dy  dx  dy  -  dxdydXdY. 

(38) 

I  I 

Assuming  z (x  ,y  )  is  a  homogeneous  random  function,  the  ensemble  average 

t  t  t  t  r  t  _  _ 

<z(x  ,y  )z(x  ,y  )>  simply  becomes  p(x,y''  where  p  is  the  surface  correlation 
function.  Equation  (38)  then  transforms  to 


I  <  E<S)H(S)*>  _  .  lli  JJJJdxdydxdY  p(x  y)  1  _  1  _ 

^  2»7  - oo  •  (X+2  x)(X"2  x) 


yj'W  4*  J.  4>  I  _  4*  4*  —  _ 

•  Z~“C  C  C  C  A  A  A  A  exp(  -  jk(  (S  +S  )  (X+0  x) - (S  +S  ) (X-0  x))) 
mnpq  mnpqmnpq  J  ''  n  m  '  2  p  q  2  ''' 
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Now  x  is  on  the  order  of  the  surface  correlation  length  (i.e.,  on  the  order  of 
centimeters).  Therefore,  it  is  negligible  relative  to  X  (which  is  many 
kilometers)  in  the  denominator  terms  of  Eq.(39).  The  imaginary  parts  of  the 
eigenangle  which  lead  to  modal  attenuation  can  also  be  ignored  when  multiplied 

by  x  in  the  exponents  of  Eq.(39).  Moreover,  the  real  parts  of  the  sines  of 
the  eigenangles  are  very  nearly  one  and  may  be  replaced  by  this  value  where 

they  multiply  x  in  Eq.(39).  Using  these  conditions  and  the  following  spectrum 
relationship,  w(k  k  ),  for  the  random  surface, 

x  >  y 

co 

2ww(kx,k  )  =  //  p(x,y)exp(-jk  x-jk  y)dxdy,  (40) 

J  —QO 


Eq.(39)  transforms  to 


\  <  E<S>H<S>*  > 


Jb_L“  ®  dXdY 

2  2nw(2k, o)//  „2 

'  -00 


zszs 

mnpq 


C  C  C*C*A  A  A*A*  exp( - jk( (S  +S  )X) - (S*+S*)X) )e  (z„)e  (zT)e*  (z„)e*  (z~)  1  . 
[mnpqmnpq  ^VJVVm  n//vp  q  ymv  T7  ynv  T7  ypv  T7  yqv  T7J 

(41) 


It  will  next  be  assumed  that  the  field  is  sufficiently  constant  over  the 
illuminated  area  so  that 


J  f (X)dX  -  AXf(X) . 


(42) 


Where  AX  is  the  pulse  length  and  X  the  range  to  the  midpoint  of  the  il¬ 
luminated  area.  The  final  integral  reduction  of  Eq.(41)  is  made  by  replacing 
dY  by  Xd<j>  and  by  assuming  an  azimuthal  beam  pattern  of  the  form 


G(<f>)  «=  exp(-in(2)  (<f>/4>Q)2) , 


(43) 


where  4>q  is  the  3dB  halfwidth  in  radians.  With  this  replacement  Eq.(41) 
reduces  to 


1<E(S)H(S)*> 

2  y  y 


f  n  nd>  w(2k  o) 

^n2^  ° 


—  XXXS 
X  mnpq 


* [CmCnCpCqAmAnApAqexp('jk((Sm+Sn)X(Sp+Sq)X))eym(zT)eyn(zT)eyp(zT)eyq(zT)] ' 


(44) 


Equation  (44)  exhibits  the  following  expected  backscatter  behavior 

1.  That  the  backscatter  power  flux  is  proportional  to  the  power 
spectrum  of  the  surface  evaluated  at  the  Bragg  wavenumber. 

2.  An  inverse  falloff  with  distance. 

3.  Proportionality  to  pulse  length  and  beamwidth. 

4.  Two  way  attenuation  between  transmitter  location  and  scatterer. 

5  Essentially  a  fourth  power  dependence  on  grazing  angle.  Recall 
that  9  is  the  angle  of  incidence  so  that  the  grazing  angle,  7, 
equals  n/2-9  and  that  to  good  approximation  Cm  -  sin7m  =*  7m. 

Though  somewhat  outdated,  the  wave  height  power  spectrum  wp(k)  suggested 
by  Pierson^  is  used  in  subsequent  calculations.  If  a  is  used  for  the  rms  sur¬ 
face  height,  w  and  Wp  are  related  as  follows 
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Mi-* 


(45) 


°2  ~  2 n  Si  w<kx-kv)dkxdkv  ”2 *  S  kdk  f  w(k,^)d^  -  J  w  (k)dk. 

-co  x  y  x  y  o  o  o  v 

If  w(k,^)  -  w(k)g(^)  with  g(^)  defined  such  that 

i  2jt 

f  gWW  -  1,  (46) 

o 


then 


00 


CO 


J  w(k)kdk  -  /  w  (k)dk, 
o  o  F 


(47) 


so  that  w(k)  -  wp(k)/k.  Equation  (44)  expressed  in  terms  of  the  Pierson 
spectrum  and  in  terms  of  grazing  angles  becomes 


<  Ey<S)Hy(s)*>  *  -  2, 


“n<2k> 


2k 


—  yyyy. 

X  mnpq 


*  [VnVqAnAmApAqexP(-jk((cos7rm+COS7n)  X  ’  (cos7*+cos7*)X) ) 


eym^zT^eyn^zT^eyp(zT^eyq^zT^] * 


(49) 


For  purposes  of  path  loss  calculations  presented  subsequently,  Equation  (49) 

2 

is  multiplied  by  the  effective  receiver  area  7r/k  and  divided  by  the  antenna 
gain,  3/2,  times  the  free  space  power  radiated  by  the  dipole.  That  is 
Equation  (49)  is  to  be  multiplied  by  the  factor. 


n  2  3  n 
k2  3  4*k 


4 


Pf. 
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Multiplying  Eq.(49)  by  this  factor,  discarding  the  negative  sign  which  simply 
signifies  energy  propagating  in  the  negative  x  direction  and  expressing  the 
result  as  path  loss  in  dB  gives 


PL(dB)  -  -10  Log10  |2.581xl06  f£h  4>Q 

*  [vnVqVmAJ  Xl  exp(-jk[(cos7m+cos7n)X  -  (cos7*+cos7*)X)]] 

*  eym(zT)  eyn(zT)  eyp(zT)eyq<zT>}-  <50> 

where  the  frequency,  f  is  to  be  expressed  in  gigahertz,  the  beamwidth  <j>Q  in 
radians,  and  all  other  quantities  in  MKS  units. 

3.0  GRAVITY  WAVE  EFFECTS 

The  gravity  wave  or  long  wavelength  region  of  the  ocean  surface  spectrum 
has  several  effects  which  are  generally  allowed  for  in  very  approximate 
fashion.  The  first  of  these  which  will  be  discussed  is  the  effect  on  the  for¬ 
ward  propagating  wave.  In  this  study  this  effect  is  approximated  by  simply 

modifying  as  follows  the  ground  reflection  coefficient  used  for  the  waveguide 
2  3 

calculations  ’ 


K<2k> 


2k 


AX  ww 
X  mnpq 


(Rim)R  "  (Rlm)F  exp(-2(kS)2  (n2(o)-/?2)),  Re(n2(o) -/?2)  >  0 
“  (Rim)F-  Re(n2(o)  -02)  <  0. 


(51) 
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Here  the  subscript  R  on  R^  signifies  the  rough  surface  value  whereas  the  sub¬ 
script  F  signifies  the  flat  earth  value.  Also,  5  is  the  rms  surface  height, 

O 

n  (o)  is  the  square  of  the  modified  refractive  index  at  the  ground  and  0  is 

2 

the  modal  propagation  constant  divided  by  wavenumber.  For  n  (o)  -  1,  the 

o  9  9 

quantity  n  (o)  -  0^  -  sin  7ffl  where  7  is  the  plane  wave  grazing  angle  at  the 
surface  and  so  the  exponent  in  the  term  in  that  limit  reduces  to  the  well 
known  factor  1/2  $  with  $  the  Rayleigh  roughness  parameter.  The  rms  surface 
height  used  in  subsequent  calculation  has  been  determined  from  the  expression 

359 

S2  =  J'w  (k)dk,  (52) 

o  ^ 

by  numerical  integration.  As  will  be  shown  later,  this  gives  values  which 

2 

agree  quite  well  with  the  often  used  formula  S~  .0051u  ,  where  u  is  windspeed 
in  m/s . 

The  composite-roughness  theory  of  sea  surface  backscatter  under  free 
space  conditions  and  within  the  line  of  sight  makes  allowance  for  grazing 
angles  of  incidence  by  introducing  a  shadowing  function,  G(y) ,  and  by  allow¬ 
ing  for  tilting  of  the  large  scale  surface.  Discussing  the  first  of  these, 
the  shadowing  concept  is  taken  over  in  the  present  study  by  simply  multiplying 
Eq.(50)  by  G(7rQ)  where  the  subscript  r  pertains  to  the  real  part  of  the  graz¬ 
ing  angle  and  the  subscript  0  pertains  to  the  least  attenuated  mode.  Using 
the  real  part  of  the  eigenangle  certainly  seems  to  be  a  reasonable  generaliza¬ 
tion  for  ducting  environments  for  which  7r»7^  (subscript  i  for  imaginary 
part).  For  standard  or  near  standard  propagation  conditions  the  assumption 
becomes  quite  suspect  since  then  7^>7r-  Nevertheless,  for  lack  of  a  rigorous 
prescription  the  assumption  has  been  used  for  all  cases  in  the  subsequent  cal¬ 
culations.  Also,  it  is  clear  that  it  would  be  preferable  to  multiply  each 
term  in  the  quadruple  sum  which  appears  in  Eq.(50)  by  an  appropriate  bistatic 
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shadowing  function.  However,  since  the  eigenspace  is  quite  densely  packed  the 
above  procedure  seems  adequate  especially  in  view  of  large  uncertainties  as¬ 
sociated  with  the  wave  height  spectrum.  The  shadowing  function  used  is  that 
9 

given  by  Wagner 

(1  +  erf(v))  (l-exo(2B)) 

G<7ro)  -  4B  *  <53) 

where 

v  -  7ro/(2f2)1/2,  (54) 


B  -  (exp(-v2)  -  Jit  verfc(v)  )/(4./irv) ,  (55) 

2  x  0 

erf(x)  -  J  exp(-t  )dt  -  1  -  erfc(x) .  (56) 

Jit  o 

Also,  f  in  Eq.  (54)  is  the  rms  slope  of  the  large  scale  surface.  Using  for- 

10  2 

mulas  given  by  Chan  and  Fung  ,  the  upwind  slope  variance,  fu,  and  crosswind 
2 

slope  variance,  fc,  are  given  by 

.359 

?u  -  5  f°k2  “p(k)<lk.  (57) 

.359 

d  -  2  J>2  “p<k>dk-  <58> 

The  upper  limit  on  the  integrals  should  be  frequency  dependent.  However,  Chan 
and  Fung  have  had  moderate  success  for  their  line  of  sight  cross  section  cal¬ 
culations  using  the  value  0.359  and  it  has  therefore  been  used  in  the  present 
calculations.  The  quantity  ft  is  given  by 


f2  .003+1. 92x10' 3  U 
c  12  5 

P  “  ~2 - Zi - L^-  (59) 

3.16x10  J  U12  5 

Where  U^2  5  is  the  wind  speed  in  meters/sec  at  an  altitude  of  12.5  m. 

Equation  (59)  was  given  by  Cox  and  Munk33  The  results  given  in  the  following 

section  presume  wind  speeds  at  19.5  m.  To  calculate  /9  for  a  given  wind  speed 

in  m/s  at  12.5  m  the  so  called  frictional  velocity,  u^,  is  first  calculated  by 

12 

iteration  from  the  equations 

ir.2„  !^f  a  (l95ol 

Oig  j  \  z  J  ’  (00) 

zQ  -  . 684/uf  +  4.28x10'5u|  -  4.43xl0'2.  (61) 

with  the  friction  velocity  determined,  U^2  ^  in  m/s  is  found  from  the  equation 

!°2  U12  5  =  It,  (i|^).  (62) 

o 

Using  this  value  in  Eq.(59)  gives  /3  which  in  turn  enables  the  calculation  of 
2  2 

fc  and  fu  as  well  as  the  corresponding  shadow  function  (53).  Results  given  in 
the  following  section  presume  upwind  conditions. 

Approximate  allowance  for  tilting  caused  by  the  large  scale  surface 
roughness  is  made  by  simply  noting  that  if  the  slope  of  the  large  scale  sur¬ 
face  is  r) ,  then  the  grazing  angle  becomes  7  +  r/ .  With  this  replacement,  each 
term  of  the  quadruple  sum  in  Eq.(50)  should  be  averaged  over  the  illuminated 

portion  of  the  surface.  However,  the  eigenangles  are  on  the  order  of  several 

f— )l/2 

milliradians  whereas  ^*7  2  J  ^  ^  is  on  the  order  of  tenths  of  radians.  Here  the 
bar  denotes  averaging  and  the  subscript  ii  denotes  that  the  averaging  is  to  be 
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( _ 7  !/2  ,  , 

taken  over  the  illuminated  portion  of  the  surface.  Because  1»?2J  ^  »  |7m|  . 

the  principal  term  in  the  quadruple  product  is  the  one  which  varies  as  Crt... 

Thus,  as  a  simplifying  assumption  each  of  the  grazing  angles  in  Eq.(50)  is  re- 

r_ a/4 

placed  by  7  + 


'ii 


This  allows  for  the  path  loss  to  be  written  as 


wp(2k)  ^ 

PL(dB)  -  -20  Log10[l.607xl03fg2  "fk“  T  G(7ro))1/2 


zs 

mn 


(7m+  )(7n+  )exp(-jk((cos7n+cos7m)X))eym(zT)eyn(zt)  j  j,  (63) 


r_iV4 


where  the  vertical  lines,  ||,  denote  absolute  values. 

Slope  averaging  over  the  illuminated  surface  has  been  discussed  by 
9  13 

Wagner  and  McDaniel  .  The  grazing  angles  dealt  with  here  are  sufficiently 
small  that  the  averaging  equation  may  be  written  as 


V  > 


U 


_2 _ 

J2n$2 


J>4ex :p(-.|2/2f2)d.>  -  3f\ 


(64) 


where  f  is  the  rms  slope  of  the  large  scale  surface  previously  discussed  in 

connection  with  Eq.(54).  For  the  purpose  of  the  calculations  given  in  the 

o 

following  section,  f  has  been  calculated  using  the  upwind  formula  given  by 
Eq . (57)  . 
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4.0  RESULTS 


Tappert  has  given  backscatter  results  for  a  single  surface  realization 
for  the  standard  atmosphere  and  for  two  evaporation  ducts.  The  latter  are 
characterized  by  duct  heights  of  14  and  28  m  and  have  been  described  by 
Hitney^  Tappert' s  calculations  are  for  a  frequency  of  9.6  GHz,  a  transmitter 
height  of  25  m,  a  pulse  width  of  0.5  ps  (in  our  convention  AX  -  cT  with  c  the 
speed  of  light  and  T  pulse  width),  and  a  3  dB  full  beam  width  of  1.2°.  In 
this  section  waveguide  backscatter  results  based  on  Eq.(63)  will  be  compared 
with  Tappert's  results.  Tables  1  and  2  below  give  the  constituents  of  Eq.(63) 
calculated  from  the  Pierson  spectrum  and  the  shadowing  function  calculated 
using  Wagner's  formula.  The  method  of  calculation  has  been  described  in  the 
previous  section. 


Table  1.  Wave  spectrum  related  quantities. 


Wind  Speed  (kts) 

JT~  (m) 

c2 

wp(2k)/2k  (m4) 

10 

.141 

.0103 

3 . 18e-4 

6 . 26e-13 

20 

.565 

.0186 

1 . 04e- 3 

1 . 63e-12 

30 

1.277 

.0306 

2 . 81e- 3 

3.09e-12 

40 

2.272 

.0444 

5 . 91e-3 

4 . 85e-12 

Table  2. 

Shadow 

function. 

Wind  Speed  (kts) 

STD  ATMOS 

14m  Duct 

28m  Duct 

10 

1. 01e-2 

6.45e-2 

9 . 91e- 2 

20 

1. 01e-2 

4 . 68e  -2 

7 . 14e  -  2 

30 

1. 01e-2 

3 . 60e-2 

5 . 45e  -  2 

40 

1. 01e-2 

2 . 95e-2 

4.47e-2 
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The  quantities  in  table  1  which  depend  upon  integration  of  the  Pierson 

spectrum  have  been  numerically  integrated  using  a  Romberg  integration  scheme 

with  a  self  adjusting  step  size.  Tabular  values  of  the  rms  bump  height  will 

2 

be  seen  to  agree  quite  well  with  the  often  used  formula  a  ~  .0051u  with  u  in 
2 

m/s.  Values  of  f  given  in  table  1  agree  quite  well  with  values  interpolated 
from  tabulations  given  by  Chan  and  Fung.  For  the  standard  atmosphere  the 
shadow  function  is  independent  of  wind  speed.  The  reason  for  this  is  that  the 
lowest  attenuated  mode  (indeed  all  modes  calculated  for  this  case)  is  evanes¬ 
cent  at  the  ground  and  so  the  rms  height  correction  to  the  ground  reflection 
coefficient  does  not  apply. 

Figures  1,  2,  and  3  show  path  loss  as  a  function  of  range  for  the  stan¬ 
dard  atmosphere,  the  14  m  duct  and  the  28  m  duct  in  that  order.  Each  plot 

shows  results  for  wind  speeds  of  10,  20,  30  and  40  knots.  The  arrow  on  the 
horizontal  axis  shows  the  normal  horizon.  Thus,  a  sufficient  number  of  modes 
have  been  used  for  the  calculations  to  allow  for  the  waveguide  fields  to  be 
carried  well  within  the  normal  horizon.  All  modes  with  attenuation  rates  less 
than  20  dB/km  were  used  for  the  standard  atmosphere  calculations  and  all  modes 
with  attenuation  rates  less  than  2.1  dB/km  were  used  for  the  evaporation  duct 
calculations.  In  particular,  24  modes  were  used  for  the  standard  atmosphere 
calculation;  between  82  and  97  modes,  depending  upon  windspeed,  were  used  for 
the  14  m  duct  calculations  and  between  94  and  109  modes  were  used  for  the  28  m 
duct  calculations.  For  ranges  beyond  and  comparable  to  the  normal  horizon 
probably  2  or  3  modes  at  the  most  would  have  been  sufficient  for  the  propaga¬ 
tion  geometry  considered.  Figures  1  through  3  show  that  the  mode  sums  break 

down  for  the  25  m  transmitter  height  at  ranges  between  about  4  and  9  km. 

Curves  for  the  standard  atmosphere,  the  14  m  duct  and  the  28  m  duct  all  show  a 
tendency  towards  saturation  as  the  wind  speed  increases.  Because  of  field 
evanescence  at  the  ground,  application  of  shadowing  concepts  is  particulary 
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*  For  the  geometry  considered,  the  extended  radio  horizon  for  the  14  m  duct  is 
about  56  km  and  about  81  km  for  the  28  m  duct. 


questionable  in  the  case  of  the  standard  atmosphere.  In  the  absence  of  al¬ 
lowance  for  shadowing,  the  path  loss  in  this  instance  would  be  reduced  by 
about  20dB  (see  table  2). 

Figures  4  through  15  show  comparisons  with  Tappert's  single  surface 
realization  calculations.  To  produce  the  comparisons,  the  path  loss  has  been 
subtracted  from  the  160  dB  system  gain  (40  dB  transmitter  gain,  40  dB  receiver 
gain  and  80  dB  transmitted  power).  In  addition  to  average  signal  levels,  the 
parabolic  equation  results  give  some  measure  of  fluctuation  about  the  average 
field.  This  presumably  results  from  the  features  of  the  illuminated  surface 
varying  with  range.  To  estimate  rms  fluctuations  with  waveguide  formalism 
would  require  carrying  the  scatter  calculation  to  higher  order.  No  attempt 
has  been  made  to  do  that.  Figures  4  through  15  show  that  the  average  field 
fall  off  (attenuation  rate)  as  determined  by  the  parabolic  equation  and 
waveguide  formalism  are  similar. 

To  give  an  idea  of  the  magnitude  of  the  difference  between  the  average 
field  as  determined  by  the  two  methods,  table  3  gives  the  difference  at  20  km 
(the  average  field  determined  by  the  parabolic  equation  calculation  has  been 
determined  by  subjectively  drawing  a  curve  through  the  center  of  the  fluctuat¬ 
ing  signal) . 


Table  3.  Comparison  of  average  fields  at  20  km. 


Wind  Speed  (kts) 


10 

20 

30 

40 


[E  -  E  ]  (dB) 
1  ywg  ype J v  ’ 


STD  Atmosphere 


14  m  Duct 


+  11.4 
+  6.0 
+  8.8 
+  11.1 


-  2.3 

-  8.0 
0.0 
0.0 


28m  Duct 

-  8.0 

-  4.6 

-  7.1 

0.0 
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The  subscripts  WG  and  PE  in  table  3  denote  waveguide  and  parabolic  equation 
results  respectively.  The  bar  over  Eypg  denotes  the  subjective  average  men¬ 
tioned  above.  The  largest  discrepancies  occur  for  the  standard  atmosphere 
with  the  waveguide  results  consistently  higher  than  the  parabolic  equation 
results.  The  most  favorable  comparison  is  for  the  14  m  duct  where  there  is 
essentially  no  difference  for  the  30  and  40  knot  wind  speed  cases.  In  con¬ 
trast  to  the  standard  atmosphere  comparisons  the  waveguide  values  are  less 
than  the  parabolic  equation  results  for  the  14m  and  28m  evaporation  ducts. 

No  attempt  has  been  made  to  use  the  same  wave  height  power  spectrum  used 
by  Tappert;  but,  rather,  in  this  initial  effort  the  decision  was  made  to  use 
the  Pierson  spectrum  which  has  been  used  by  Chan  and  Fung,  and  McDaniel  with 
modest  success  at  more  oblique  incidence.  Of  lesser  importance  is  the  fact 
that  Tappert 's  results  are  for  an  infinitely  conducting  ground,  whereas  the 
present  calculations  are  for  a  seawater  conductivity  of  21  Si/m  and  a  seawater 
dielectric  constant  of  54.3.  In  view  of  the  approximate  allowance  for  surface 
roughness  on  forward  propagation,  the  non-rigorous  introduction  of  shadowing 
and  tilting  concepts  along  with  the  omission  of  multiple  scatter  effects,  the 
agreement  exhibited  by  figures  4  through  15  and  table  3  is  very  much  better 
than  expected. 


5.0  DISCUSSION 

Generally,  normal  mode  methods  have  been  used  to  assess  the  adequacy  of 
parabolic  equation  solvers  for  propagation  in  either  acoustic  or  electromag¬ 
netic  channels.  Recent  results  by  Tappert,  using  parabolic  equation  and  Monte 
Carlo  methods,  for  electromagnetic  backscatter  in  evaporation  ducts  afford  an 
opportunity  to  use  parabolic  equation  results  to  assess  the  adequacy  of  using 
waveguide  methods  for  treating  backscatter  for  very  grazing  angles  (i.e.  modal 
angles  on  the  order  of  several  milliradians) .  In  particular,  the  parabolic 
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equation  results  allow  for  an  assessment  of  classical  approximations  used  to 
treat  long  surface  wavelength  effects.  These  include  applicability  of 
Kirchhoff  theory  for  forward  propagation  along  with  shadowing  and  tilting  of 
the  local  surface.  In  this  study,  results  at  9.6  GHz  have  been  generated 
based  on  first  order  Bragg  scatter  in  a  waveguide  environment  coupled  with  the 
long  surface  wavelength  effects  just  mentioned.  Comparisons  of  the  results 
with  Tappert's  single  surface  realization  are  shown  in  figures  4  through  15. 
Those  comparisons  show  very  similar  attenuation  rate  behavior  for  the  ducting 
environments  considered  (i.e.  normal  atmosphere,  14  and  28  m  evaporation  ducts 
with  10,  20,  30,  and  40  knot  wind  speeds).  Average  levels  compare  to  within 
about  ±  1.0  dB  (see  table  3)  .  Some  of  this  difference  is  probably  due  to  the 
fact  that  Pierson's  surface  height  spectrum  has  been  used  in  the  present  study 
whereas  Tappert  used  yet  another  spectrum.  Nevertheless,  in  view  of  the  crude 
and  certainly  questionable  methods  for  treating  long  wavelength  effects,  the 
comparisons  are  surprisingly  good  and,  at  the  very  least,  intriguingly 
curious.  Certainly,  additional  comparisons  should  be  made. 

The  object  of  this  study  has  been  to  explore  the  possibility  of  using 
waveguide  concepts  for  backscatter  calculations  and  not  the  development  of  a 
general  purpose  computer  program.  Thus,  surface  quantities  such  as  rms  bump 
height,  tilting  angle,  shadowing  function,  etc.  have  been  calculated  using 
routines  separate  from  the  waveguide  program.  An  automated  program  requires 
the  merging  of  these  programs.  Formulas  for  line  of  sight  backscatter  for 
horizontally  stratified  environments  should  be  developed  and  added  to  the 
program.  This  would  allow  scatter  calculations  for  ranges  less  than  »  8  km 
where  the  mode  sums  break  down  in  the  present  case  studies.  Only  horizontal 
polarization  has  been  treated.  Expectations  are  that  vertical  polarization 
would  behave  in  similar  fashion  beyond  the  horizon.  Within  the  line  of  sight, 
however,  differences  are  likely  and  results  applicable  to  vertical  polariza¬ 
tion  should  be  developed  lastly,  first  order  scatter  results  were  used  in 
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the  present  study.  To  estimate  variances  of  the  scattered  fields  higher  order 
results  are  required.  That  would  be  a  challenging  and  interesting  extension 
of  the  present  study. 
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RANGE  (KM) 

FREQUENCY**  9600.0  MHZ  Figure  1.  Path  loss  vs.  range  for  the  standard  atmosphere. 

Transmitter-recei ver  height  25m. 
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Figure  2.  Path  loss  vs.  range  for  a  14m  evaporation  duct. 
Transmitter-receiver  height  25m. 
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Figure  4.  Clutter  power  vs.  range.  Comparison  between  pe  and  waveguide 
results.  Standard  atmosphere,  10  knot  wind. 
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Clutter  power  vs.  range.  Comparison  between  pe  and  waveguide 
results.  Standard  atmosphere,  20  knot  wind. 
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Figure  6.  Clutter  power  vs.  range.  Comparison  between  pe  and  waveguide 
results.  Standard  atmosphere,  30  knot  wind. 
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Figure  7.  Clutter  power  vs.  range.  Comparison  between  pe  and  waveguide 
results.  Standard  atmosphere,  40  knot  wind. 
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Range  (km) 


Figure  8.  Clutter  power  vs.  range.  Comparison  between  pe  and  waveguide 
results.  Fourteen  meter  evaporation  duct,  10  knot  wind. 
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Figure  9.  Clutter  power  vs.  range.  Comparison  between  pe  and  waveguide 
results.  Fourteen  meter  evaporation  duct,  20  knot  wind. 
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Figure  11.  Clutter  power  vs.  range.  Comparison  between  pe  and  waveguide 
results.  Fourteen  meter  evaporation  duct,  40  knot  wind. 


Figure  12.  Clutter  power  vs.  range.  Comparison  between 
results.  Twenty  eight  meter  evaporation  due 
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Figure  13.  Clutter  power  vs.  range.  Comparison  between  pe  and  waveguide 
results.  Twenty  eight  meter  evaporation  duct,  20  knot  wind. 
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Figure  14.  Clutter  power  vs.  range.  Comparison  between  pe  and  waveguide 
results.  Twenty  eight  meter  evaporation  duct,  30  knot  wind. 
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